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Abstract 

It has been long-conjectured that the crossing number cv{K m , n ) of the complete bi- 
partite graph K m , n equals the Zarankiewicz Number Z(m,n) := L 22 ^] LifJ \ !1 ~T~\ Lf J- 
Another long-standing conjecture states that the crossing number cr(_ftT„) of the com- 
plete graph K n equals Z(n) := \ \\ \ [^^-J \ L 2 ^] ■ 111 this paper we show the 
following improved bounds on the asymptotic ratios of these crossing numbers and 
their conjectured values: 

(i) for each fixed m > 9, linin^oo cr(K m ^ n )/Z(m, n) > 0.83m/(m— 1); 

(ii) Um n _Kx, cr(K n n )/Z(n, n) > 0.83; and 
(hi) lim n ^oc cx{K n )/Z{n) > 0.83. 

The previous best known lower bounds were 0.8m/ (m — 1),0.8, and 0.8, respectively. 
These improved bounds are obtained as a consequence of the new bound cr(K 7t7l ) > 
2.1796n 2 - 4.5n. To obtain this improved lower bound for cr(7^7 ira ), we use some el- 
ementary topological facts on drawings of K-zj to set up a quadratic program on 6! 
variables whose minimum p satisfies cr(.KY )n ) > (p/2)n 2 — 4.5n, and then use state-of- 
the-art quadratic optimization techniques combined with a bit of invariant theory of 
permutation groups to show that p > 4.3593. 
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1 Introduction 



In the earliest known instance of a crossing number question, Paul Turan raised the prob- 
lem of calculating the crossing number of the complete bipartite graphs K miTl . Turan's 
interesting account of the origin of this problem can be found in 



We recall that in a drawing of a graph in the plane, different vertices are drawn as different 
points, and each edge is drawn as a simple arc whose endpoints coincide with the drawings 
of the endvertices of the edge. Furthermore, the interior of the arc for an edge is disjoint 
from all the vertex points. We often make no distinction between a graph object, such as 
a vertex, edge, or cycle, and the subset of the plane that represents it in a drawing of the 
graph. 

The crossing number cr(G) of a graph G is the minimum number of pairwise intersections 
of edges (at a point other than a vertex) in a drawing of G in the plane. 

Exact crossing numbers of graphs are in general very difficult to compute. Long-standing 
conjectures involve the crossing numbers of interesting families of graphs, such as Km,n 
and K n . On a positive note, it was recently proved by Glebskii and Salazar [Hj that the 
crossing number of the Cartesian product C m x C n of the cycles of sizes m and n equals its 
long-conjectured value, namely (m — 2)n, at least for n > m{m + 1). For a recent survey 
of crossing number results, see for instance Shahrokhi et al. j2Uj . 

= Z(m, n) for all 
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Zarankiewicz published a paper |25j . in which he claimed that cr(K„ 
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However, several years later Ringel and Kainen independently found a hiatus in Zaran- 
kiewicz's argument. A comprehensive account of the history of the problem, including a 
discussion of the gap in Zarankiewicz's argument, is given by Guy 10 . 

Figure Q shows a drawing of with 8 crossings. As Zarankiewicz observed, such a 
drawing strategy can be naturally generalized to construct, for any positive integers m,n, 
drawings of K m ^ n with exactly Z(m, n) crossings. This observation implies the following 
well-known upper bound for cr(K mtn ): 

cr(K mi „) < Z(m,n). 

No one has yet exhibited a drawing of any K mn with fewer than Z(m, n) crossings. In 
allusion to Zarankiewicz's failed attempt to prove that this is the crossing number of K m ^ n , 
the following is commonly known as Zarankiewicz's Crossing-Number Conjecture: 



cx(K v 
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for all positive integers m,n. 
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Figure 1: A drawing of with 8 crossings. A similar strategy can be used to construct 
drawings of K m>n with exactly Z(m,n) crossings. 



In 1973, Guy and Erdos wrote "[Al]most all questions that one can ask about crossing 
numbers remain unsolved" More than three decades later, despite some definite progress 
in our understanding of this elusive parameter, most of the fundamental and more important 
questions about crossing numbers remain open. Zarankiewicz's Conjecture has been verified 
by Kleitman jllj for min{m, n} < 6 and by Woodall 24 for the special cases 7 < m < 8, 
7 < n < 10. 

Since the crossing number of K m ^ n is unknown for all other values of m and n, it is natural 
to ask what are the best general lower bounds known for cr^K mn ). A standard counting 
argument, together with the fact that cx{K^ jn ) is as conjectured, yields the best general 
lower bound (j2J) known for cx(K m n ). It goes as follows: suppose we know a lower bound 
c r on cx(K r>n ), for 2 < r < m < n. Each crossing in the embedding of K m ,n lies m C^lT^) 
distinct K r ^ n C K m ^ n . As there are in total (™) distinct K r>n 's, one obtains 

c,J m ) 

cr (K m;n ) > , for r = 5 one derives cr(if m>n ) > 0.8 Z(m, n). (2) 

( r-2 ) 

A small improvement on the 0.8 factor (roughly to something around 0.8001) was recently 
reported by Nahas |15j . 

Zarankiewicz's Conjecture for Kj )n states that 

( 2.25n 2 -4.5n + 2.25, n odd, n>7 
1 2.25n 2 — 4.5n, n even, n > 8 

As we observed above, this has been verified only for n = 7, 8, 9, and 10. Using cr(i^7 io) = 
180, a standard counting argument gives the best known lower bounds for cr(Kj tn ) for 
11 < n < 22. However, for n > 23, the best known lower bounds for cr(K-/ tn ) are obtained 
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by the same counting argument, but using the known value of ci{K^^ n ) instead of cr(i£Y 10). 
Summarizing, previous to this paper, the best known lower bounds for cr(7, n) were: 



An elementary calculation shows that this is an improvement, for all n > 23, on the bounds 
for cr(Kj tn ) given in (J3J). 

The strategy of the proof can be briefly outlined as follows. Let (A, B) be the bipartition 
of the vertex set of K^ >n , where \A\ = 7 and \B\ = n > 2. Let b, b' be vertices in B. In any 
drawing T> of K^ tTl , the number of crossings that involve an edge incident with b and an edge 
incident with b' is bounded from below by a function of the cyclic rotation schemes of b and 
b' . This elementary topological observation on drawings of K2 7 naturally yields a standard 
quadratic (minimization) program whose minimum p satisfies cr(i^7 jn ) > (p/2)n 2 — 4.5n 
(see Lemma OJ)- We then use state-of-the-art quadratic programming techniques to show 
that p > 4.3593 (see Lemma 01, thus implying Theorem Q 

The rest of this paper is organized as follows. In Section 01 we review some elementary 
topological observations about drawings of K2 >n , and use these facts to set up the quadratic 
program mentioned in the previous paragraph. The bound for ci(Kj tn ) in terms of the 
minimum of this quadratic program is the content of Lemma 01 In Section 01 we prove 
Proposition 01 which gives a lower bound for the quadratic program. As we observe at the 
end of Section 01 Theorem ^ is an obvious consequence of Lemmas 01 and 01 In Section 21 
we discuss consequences of Theorem ^ the improved bound for cr(Kj tn ) implies improved 
asymptotic bounds for the crossing numbers of cx(K mn ) and cr(K n ). 



2 Quadratic optimization problem yielding a lower bound 

for cr(K m>n ) 



Our goal in this section is to establish Lemma 01 a statement that gives a lower bound 
for cv{K m ^ n ), for m < n, (and thus for cr(K^ tn )) in terms of the solution of a quadratic 
minimization problem on (m — 1)! variables. 

Let n > m be fixed. Let V denote the vertex set of if m , n , and let (A, B) denote the 




11 < n < 22 
odd n > 23 
even n > 24 



(3) 



In this paper we prove the following. 



Theorem 1. For all integers n 



cr(K 7 , n ) > 2.1796n 2 - 4.5n. 
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Figure 2: Here m = 7. Vertices bi and 6j have cyclic orderings (0134526) and (0265341), 
respectively (we write i for a% for the sake of brevity). It is easy to check that the mini- 
mum number of interchanges among adjacent elements in (0134526) required to obtain 
(0265341)- 1 (namely (0143562)) is 2. Thus, Q((0 1 3 4 5 2 6), (0 1 3 4 5 2 6)) = 2. There- 
fore, there must be at least two crossings (as it is indeed the case in the drawing above) 
that involve edges incident with bi and by 



bipartition of V such that each vertex of A = {ao, a\, . . . , a m _i} is adjacent to each vertex 
of B = {6 Q) 6i, . . . ,6 n _i}. 

Consider a fixed drawing T> of K m>n . To each vertex bi we associate a cyclic ordering 7Tx>(&j) 
of the elements in A, defined by the (clockwise) cyclic order in which the edges incident 
with bi leave 6j towards the vertices in A (see Figure Let II denote the set of all cyclic 
orderings of {ao, «i, . . . , a m _i}. Note that | XT | = m\/m = (m — 1)1. 

Following Kleitman let crx>(bi,bj) denote the number of crossings in T> that involve an 
edge incident with bi and an edge incident with bj. Further, let p\ y pi € II and Q(pi,p2) be 
the minimum number of interchanges of adjacent elements of p\ required to produce p^ . 
Then, for all bi, bj with 6, ^ bj, 

ctvih, bj) > Q(Trv(bi),7rv(bj)). (4) 

This inequality is stated in jllj and proved in |24j . This observation alone yields a lower 
bound for cr(K m n ), as follows. Fix any drawing T> of K m ^ n . For each p £ II, let 

1 

x p := - \{bi £ B | nv{bi) = p}\ 
n 

The matrix Q can be viewed as the matrix of quadratic form Q(-, ■) on the space lRl n L It 
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follows from (JIJ) that 



cr(V)> Q(p,p')(x P n)(x p/ n) + J2Q(PiP)( % 
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where we have used the (easily verifiable, see e.g. |24p fact that Q(p, p) = \m/2\ |_(m — 1) /2j 
for every p G II. 

Since the drawing T> was arbitrary, we have proved the following. 

Lemma 2. Let Q be the (m — 1)1 x (m — 1)! matrix of the form Q(-, •), and let e denote the 
all ones vector. Then, for every integer n > m > 2, 



cr(K m>n ) > — ( nmm{x T Qx \ x G R^™ ,e T x = 1} 



m 




m — 1 






2 



cr(K-j^ n ) > — ( nmin^ x T Qx 



x eRf,e T x = l} - 9) . 



Remark For obvious reasons {Q is a 720 x 720-matrix) we do not include in this paper the 
matrix Q in table form. As we mentioned above, Q(p, p) = 9 for every p £ II, and therefore 
all the diagonal entries of Q are 9. It is not difficult to show that Q(p,p') < 8 if p ^ p' , so 
every non-diagonal entry of Q is at most 8. The calculation of the entries of Q, using the 
definition of Q(-, •), and taking its symmetries into account (see Section [3 .21 on this) takes 
only a few seconds of computer time. 



3 Finding a lower bound for the optimization problem 

Our aim in this section is to find a (reasonably good) lower bound for the quadratic pro- 
gramming problem with m = 7 given in Lemma in order to obtain a (reasonably good) 
lower bound for cr(Ki tn ). The main result in this section is the following. 

Proposition 3. Let Q be the 6! x 6! matrix of the quadratic form Q(-, ■). Then 



min< x Qx 



x <E M. 6 l,e T x = l[ > 4.3593. 



We devote this section to the proof of Proposition |21 It involves computer calculations; 
more details on this are given in Section f3.8l 
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3.1 The standard quadratic programming problem 



The problem we have formulated is known as standard quadratic optimization problem. The 
standard quadratic optimization problem (standard QP) is to find the global minimizers 
of a quadratic form over the standard simplex, i.e. we consider the global optimization 
problem 

p := mmx T Qx (5) 
— xgA 

where Q is an arbitrary symmetric d x d matrix, e is the all ones vector, and A is the 
standard simplex in R d , 

A = {x € K+ : e T x = 1}. 

We will now reformulate the standard QP problem as a convex optimization problem in 
conic form. First, we will review the relevant convex cones as well as the duality theory of 
conic optimization. We define the following convex cones. 

• The d x d symmetric matrices: 
S d = {X € R d x R d , X = X T }; 

• The d x d symmetric positive semidefinite matrices: 
S+ = {XG S d , y T Xy > Vy € R d }; 

• The d x d symmetric copositive matrices: 
C d ={X€ S d , y T Xy > Vy e R d , y > 0}; 

• The d x d symmetric completely positive matrices: 

C d = {x = ^ =1 y iy f, y t GR d , Vi > (i = 1, . . . , k)}; 

• The d x d symmetric nonnegative matrices: 
M d = {X G S d , Xij > (i, j = 1, . . . , d)}. 



Recall that the completely positive cone is the dual of the copositive cone, and that the 
nonnegative and semidefinite cones are self-dual for the inner product (X, Y) := Tr(XY), 
where 'Tr' denotes the trace operator. 

For a given cone K, d and its dual cone IC* d we define the primal and dual pair of conic linear 
programs: 

p* := inf {Tr(CX) | Tr(^X) = h (i = 1, . . . , M)} (P) 

AG/C d 

d*:= sup \b T y\ y^ yi A l + S = C, S € KrA . (D) 

y^ Rm { i=i J 

If K d = St we refer to semidefinite programming, if K, d = M d to linear programming, and 
if K. d = Cd to copositive programming. 
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The well-known conic duality theorem (see e.g. Renegar ^7j) gives the duality relations 
between (P) and (D). 

Theorem 4 (Conic duality theorem). If there exists an interior feasible solution X° E 
int{K,(i) of (P), and a feasible solution of (D) then p* = d* and the supremum in (D) is 
attained. Similarly, if there exist feasible y°,S° for (D) where S° E int(K,V), and a feasible 
solution of (P), then p* = d* and the infimum in (P) is attained. 

Optimization over the cones <Sj~ and Md can be done in polynomial time (to compute an 
e-optimal solution), but some NP-hard problems can be formulated as copositive programs, 
see e.g. De Klerk and Pasechnik |12| . 

3.1.1 Convex reformulation of the standard QP 

We rewrite problem © in the following way: 

p := minTr((5xx T ). 
— xeA 

Now we define the cone of matrices 

K = {X E S d : X = xx T , x > 0} . 

Note that the requirement x E A corresponds to X E K, with Tr (ee r X) = 1. 
We arrive at the following reformulation of problem 

p = min {Tt(QX) : Tr (ee T X) = 1, X E K) . (6) 

The last step is to replace the cone tC by its convex hull, which is simply the cone of 
completely positive matrices, i.e. 

conv(£)=C2 = j* = X>yf, yteW 1 , yi>0(i = l,...,A;)|. 

Replacing the feasible set by its convex hull does not change the optimal value of problem © , 
since its objective function is linear. Thus we obtain the well-known convex reformulation: 

p = min (Tr(QX) | Tr (ee T X) = 1, X E C* d ] . (7) 
The dual problem takes the form: 

p = max {t \ Q — tee T E Q} , (8) 

where Cd is the cone of copositive matrices, as before. Note that both problems have the 
same optimal value, in view of the conic duality theorem. 
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3.2 Exploiting group symmetries 



We can reduce the number of variables in the optimization problems in (J7I8|) considerably 
by exploiting the invariance properties of the quadratic function x T Qx. This will also prove 
to be computationally necessary for the problems we intend to solve. 

Consider the situation where the matrix Q is invariant under the action of a group G of 
order k = \G\ of permutation matrices P G G, in the sense that 



Q = P T QP for all P G G. 



Then we have 
V 



min {Tr{QX) | Tr (ee T X) = 1, X G C* d ) 
min {Tr (P T QPX) \ Tr (Pee T PX) = 1, X G C* d ) for any P G G 
min {Tr (QP T XP) | Tr (ee T P T XP) =1, X G C*j} for any P G G 
1 



min < Tr I Q 



k 



Y,p t xp 



PeG 



) I Tr ( egT 1 E pTxp J = i, * e c^J 



We can therefore restrict the optimization to the subset of the feasible set obtained by 
replacing each feasible X by the group average ^ X^PeG P T XP, i.e. replacing X by its 
image under what is known in invariant theory as the Reynolds operator. Note that if 
X G C* d , then so is its image under the group average. 

In particular, we wish to compute a basis for the so-called fixed point sub space: 

A:=Iy eS d \Y = ±Yl pTxP > X e S d \ ■ 
{ PeG J 

Note that Q and ee T are elements of A (set X = Q, respectively, X = ee T ). Hence 
Q — tee T G A for any t, and 

p = max {t \ Q — tee T G C d } = max {t \ Q — tee T G^fl .4} . 

The right-hand side here is the dual of the primal problem when it was restricted to A as 
above. 

The next step is to compute a basis for the subspace A. 



3.3 Computing a basis for the fixed point subspace 

We assume for simplicity that G acts transitively as a permutation group on the standard 
basis vectors. (This holds in our setting. A more general, and computationally less efficient, 
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setting, can be found in Gatermann and Parrilo [Sj.) The theory here is well-known, and 
goes back to Burnside, Schur and Wielandt. See e.g. Cameron [H] for details. Although we 
need a basis of A, the subspace of symmetric matrices fixed by G, it is more natural to 
compute the basis X of the subspace B of all fixed by G matrices, and then pass on to A. 

The dimension of B equals the number r of orbits of G on the Cartesian square of the 
standard basis. The set of the latter orbits, also known as 2-orbits, naturally corresponds 
to certain set X of d x d zero-one matrices. Namely, for each X G X one has X^ = 1 if 
and only if Xp^^p^ = 1 for all P G G and all 1 < i < j < \U\. As G is transitive on the 
standard basis vectors, the identity matrix I belongs to X. We also have ^2 x &x = ee - T ■ 

As X is closed under the matrix transposition, i.e. X T G X for any X G X, 



Since Q G A, we will write Q = X^i=i foAi. 

It is worth mentioning that algebraically the vector space B behaves very nicely: it is closed 
under multiplication. In other words, B is a matrix algebra of dimension r, also known as 
the centralizer ring of the permutation group G. 

We proceed to describe G and B in our case. For us G is isomorphic to the direct product 
Sym(m) x Sym(2), of symmetric groups Sym(m) and Sym{2), where Sym(m) acts (as a 
permutation group) by conjugation on the d = (m — 1)! elements of II, and Sym{2) acts 
(as a permutation group) on II by switching it G II with tt~ 1 G II. 

Computing X is an elementary combinatorial procedure, that can be found in one form or 
another in many computer algebra systems, so one does not have to program this again. 
First, the permutations that generate Sym{m) x Sym{2) in its action on II are computed. 
The action of Sym{2) is already known, and is described by the permutation go, say. In its 
usual action on m symbols, Sym(m) is generated by hi = (0, 1, . . . , m — 1) and h% = (0, 1). 
These hi (for i = 1,2) act on II by mapping each tt G II to hiirh^ 1 . Denote by gi (for 
i = 1, 2) the permutations of II that realize these actions. 

Next, one computes the orbits of the permutation group Sym(m) x Sym(2) = {go, (71,52) on 
the Cartesian square II x II of II, by "spinning" (7Tj, ttj) G II x II: begin with Sij = {(vrj, ttj)} 
and apply the generators gi, < i < 2, in a loop until Sij stops growing. Then one sets 
II := II — Sij and repeats until II is exhausted. 

When m = 7 one has r = 78 and M = 56. Note that here the algebra B is not commutative. 
When m = 5 one has r = M = 6, and B commutative. 



X A = {A x ,..., A M ] = {X I X = X T G X}U{X + X T \ X G X, X ^ X T } 

is a basis of A. Each A G Xj± is a symmetric zero-one matrix and Xmsa^ ^ 
Moreover 




10 



3.4 Reformulation of the optimization problem 



We can now reformulate the dual problem by using the basis of A to obtain: 

p = max {t | Q - tee T eCdHA} = max < t \ ^J(6* - t)Ai G C d > . 

We will now proceed to derive a lower bound on p by solving the dual problem approxi- 
mately. 



3.5 Approximations of the copositive cone 

The problem of determining whether a matrix is not copositive is NP-complete, as shown 
by Murty and Kabadi ^1]. We therefore wish to replace the copositive cone Cd by a conic 
subset, in such a way that the resulting optimization problem becomes tractable. We can 
represent the copositivity requirement for a (d x d) symmetric matrix S as 

d 

P(x) := (x o x) T S(x o x) = Sijxjx] > 0, Vx € R d , (9) 

where 'o' indicates the componentwise (Hadamard) product. We therefore wish to know 
whether the polynomial P(x) is nonnegative for all x £ R . Although one apparently cannot 
answer this question in polynomial time in general, as it is an NP-hard problem, one can 
decide using semidefinite programming whether P(x) can be written as a sum of squares. 

Parrilo ^B] showed that P(x) in @ allows a sum of squares decomposition if and only if 
S G + Afd, which is a well-known sufficient condition for copositivity. Set /C|] to be the 
convex cone K? d = <Sj" + J\fd- 



Higher order sufficient conditions can be derived by considering the polynomial: 

PW(s) = P(x) (y>A = ( V S ;r rj.rj] (y>A , (10) 




and asking whether P^\x) - which is a homogeneous polynomial of degree 2(1 + 2) - has 
a sum of squares decomposition, or whether it only has nonnegative coefficients. 

For £ = 1, Parrilo ^H] showed that a sum of squares decomposition exists if and only if 1 
1 ln fact, Parrilo 1101 only proved the 'if'-part; the converse is proven in Bomze and De Klerk 4 . 
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the following system of linear matrix inequalities has a solution: 

S-S® G «S+, i = l,...,d (11) 

S$ = 0, i = l,...,d (12) 

sV+2S$ = 0, i + j (13) 

> 0, i<j<fc, (14) 

where (i = 1, . . . , d) are symmetric matrices. Similar to the I = case, we define lC d as 
the (convex) cone of matrices S for which the above system has a solution. 



We will consider the lower bounds we get by replacing the copositive cone by either JC d or 
P >Pi ■= max|t | Q -tee T G /C^} , £ G {0, 1}. (15) 



3.6 Approximations (relaxations) of the copositive cone 

We will now study the relaxation obtained by replacing the copositive cone by its proper 
subset K. d . In other words, we study the relaxation: 

= maxjt | J2(bi -t)Ai G C d | 
> po :=max|t | ^(fej-t^ e JC° d = S+ + Afjj 
We rewrite £^(6* ~ t) A i G ^° as 

M MM M M 

- t)Ai = ^2 yiAi + ^2 z i A ii where ^ j/jAj G <Sj~ and ^ ^,4; G M d - 

i=l i=l i=l i=l i=l 

Note that, since the AiS are 0-1 matrices that sum to ee T , it follows that z,- L > 0. Moreover, 

h — t = yi + Zi implies bi — t — yi > 0. 
We obtain the relaxation: 



p = max|t \ bi-t-yi>0{i=l,...,M), ^y^G<S+|. (16) 
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3.7 Block factorization 



The next step in reducing the problem size is to perform a similarity transformation that 
simultaneously block-diagonalizes the matrices A\, . . . , Am- In particular, we want to find 
an orthogonal matrix V such that the matrices 

Af-VAtV- 1 i = l,...,M, 

all have the same block diagonal structure, and the maximum block size is as small as 
possible. Note that the conjugation preserves spectra, and orthogonality of V preserves 
symmetry. 

This will further reduce the size of the relaxation Q16[) via 



Po = max 



max 



max 



^t\bi-t-y t >0(i = l,...,M),J2 Vi A i G S d j 

ji \ b i -t-y i >0(i = l,...,M), YjjiVAiV" 1 eS+j 

|i I 6 i -t-y i >0(t = l,...,M), ^?/iAiG«5+|. 



The necessity to restrict to orthogonal V's lies in the fact that there is currently no software 
(or algorithms) available that would be able to deal with non-symmetric Aj's. 

Computing the finest possible block decomposition (this would mean finding explicitly the 
orthogonal bases for the irreducible submodules of the natural module of G in its action by 
the matrices P) is computationally not easy, especially due to the orthogonality requirement 
on V . We restricted ourselves to decomposing into two blocks of equal size | x 5. Namely, 
each row corresponds to a cyclic permutation g G II, and the natural pairing {g,g~ l ) can 
be used to construct V = ^V, as follows: 

• the first half of the rows of V' are characteristic vectors of the 2-subsets {flsg -1 }, 

9 en ; 

• the second half of the rows of V 1 consists of "twisted" rows from the first half: namely 
one of the two Is is replaced by -1. 

It is obvious that V'V' T = 21 and thus V is orthogonal. 



Remark It is worth mentioning that in ^S] Schrijver essentially dealt, in a different 
context, with a similar setup, except that in his case the elements of the basis X of B were 
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symmetric and (hence) the algebra B commutative. In such a situation the elements of X 
can be simultaneously diagonalized, and the corresponding optimization problem becomes 
a linear programming problem. 

3.8 Computational results: proof of Theorem OH 

The combinatorial/group theoretic part of the computations, namely of the Aj's, V, and 
Q = was performed using a computer algebra system GAP [7] version 4.3, and its 

shared package GRAPE by Soicher [2TJ . Semidefinite programs were solved using SeDuMi 
by Sturm version 1.05 under Matlab 6.5. The biggest SDP took about 10 minutes of 
CPU time of a Pentium 4 with 1 GB of RAM. 

In addition, the results were verified using Maple. Namely, for t = po and y, the variables 
computed upon solving (j!6ft . we checked that the corresponding (matrix and scalar) in- 
equalities in (fH)|) hold. As po is a lower bound on p, we thus validated the computed value 
of po independently of the SDP solver used. 

For the test case of K^^ n we solved the relaxed problem (j!5|) with £ = 1 to obtain 

Pl « 1.9544, that is cr(K 5 , n ) > i(1.9544)n 2 = 0.9772n 2 , 

asymptotically. The correct asymptotic value is known to be ci(K^ n ) = n 2 , which shows 
the quality of the bound. The weaker bound for i = in ()15|) yields, still quite tight 

Pl pa 1.94721, that is cr(iT 5)n ) > ^(1.94721)n 2 = 0.973605n 2 . 

For the case K? }Tl we solved the relaxed problem (|15|) with £ = to obtain 

Po ~ 4.3593, that is cr^J > ^(4.3593)n 2 = 2.1796n 2 , 

asymptotically. 

Proof of Theorem ^ For the sake of completeness, we close this section with the ob- 
servation that Theorem n has been proved. It follows from Lemmas 121 and 01 □ 

4 Improved bounds for the crossing numbers of K m ^ n and K n 

Perhaps the most appealing consequence of our improved bound for cr(Ki tn ) is that it also 
allows us to give improved lower bounds for the crossing numbers of K m ^ n and K n . The 
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quality of the new bounds is perhaps best appreciated in terms of the following asymptotic 
parameters: 



A(m) 



lim 



cr(if m , n ) 



B := 



cx{K n . 



lim 

n— >oo Z(n,n) 



n-*oo Z(m, n) ' 

(see Richter and Thomassen [18 ). These natural parameters give us a good idea of our 
current standing with respect to Zarankiewicz's Conjecture. It is not difficult to show that 
A{m) (for every integer m > 3) and B both exist |18j . 

Previous to the new bound we report in Theorem^ the best known lower bounds for A(m) 
and B were A(m) > 0.8^^- and (consequently) B > 0.8. Both bounds were obtained by 
using the known value of cr(K^ tU ), and applying a standard counting argument. 

By applying the same counting argument, but instead using the bound given by Theorem^ 
we improve these asymptotic quotients to A(m) > 0.83^^- and B > 0.83. 

The improved lower bound for B has an additional, important application. It has been long 
conjectured that ci(K n ) = Z(n), where 



Z(n) = - 
V 1 4 



n 



n 



n 



but this has been verified only for n < 10 (see for instance jHj). As we did with K m n , it is 
natural to inquire about the asymptotic parameter 



C 



lim 

n^oo 



Z{n) 



In |18j it is proved that C exists, and, moreover, C > B. In view of this, our improved 
lower bound for B yields C > 0.83. 

We summarize these results in the following statement. 
Theorem 5. With Z(m,n) and Z(n) as above, 

lim Cr{Km > n } > 0.83^, lim ^4 > 0-83, and lim ^ > 0.83 □ 



We close this section with a few words on some important recent developments involving 
the rectilinear crossing number of K n . 

The rectilinear crossing number cr(G) of a graph G is the minimum number of pairwise 
intersections of edges in a drawing of G in the plane, with the additional restriction that 
all edges of G must be drawn as straight segments. 

It is known that cr(i^ n ) and cr(K n ) may be different (for instance, cr(K$) = 19, whereas 
cr(i"Cg) = 18). While we have a (non-rectilinear) way of drawing K n that shows cr(K n ) < 
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Z(n) (equality is conjectured to hold, as we observed above), good upper bounds for cx(K n ) 
are notoriously difficult to obtain. Currently, the best upper bound known is cr(K n ) < 
0.3807(2) ( see Aichholzer et al. 0). 

For many years the best lower bounds known for cr(K n ) were considerably smaller (around 
0.32 Q)) than the best upper bounds available (currently around 0.380(2))- However, re- 
markably better lower bounds have been recently proved independently by Abrego and 
Fernandez-Merchant [Tj and by Lovasz et al. ^3], and refined by Balogh and Salazar |3J- 
In pQ, the technique of allowable sequences is used to show that cf(K n ) > 0.375(2)- Lovasz 
et al. use similar methods to prove cr(K n ) > 0.37501(2) + 0(n 3 ). Recently, Balogh and 
Salazar improved this to cr(K n ) > 0.37553(2) 0(n 3 ) [3J. The importance of establish- 
ing that cr(K n ) is strictly greater than 0.375(2) + 0(n 3 ) is that it effectively shows that 
the ordinary and the rectilinear crossing numbers of K n are different in the asymptotically 
relevant term, namely n 4 . 

Acknowledgements. Etienne de Klerk would like to thank Pablo Parrilo for his valuable 
comments. 
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